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ABSTRACT 

We present new relationships between halo masses {Mi,) and several galaxy properties, including r*-band 
luminosities (L r ), stellar (M stal ) and baryonic masses, stellar velocity dispersions (er), and black hole masses 
C^bh)- Approximate analytic expressions are given. In the galaxy halo mass range 3 x 10 10 M Q < M/, < 
3 x 1O 13 M the Mh-L r , M slM -Mh, and Mbh-Mi, are well represented by a double power law, with a break 
at M/,,break ~ 3 x 10 n M Q , corresponding to a mass in stars M star ~ 1.2 x 10 10 M©, to a r*-band luminosity 
L, - ~ 5 x 10 9 L Q , to a stellar velocity dispersion a ~ 88 km s" 1 , and to a black hole massM^// ~9x 1O 6 M . The 
<7— Mi, relation can be approximated by a single power law, though a double power law is a better representation. 
Although there are significant systematic errors associated to our method, the derived relationships are in good 
agreement with the available observational data and have comparable uncertainties. We interpret these relations 
in terms of the effect of feedback from supernovae and from the active nucleus on the interstellar medium. We 
argue that the break of the power laws occurs at a mass which marks the transition between the dominance of 
the stellar and the AGN feedback. 

Subject headings: cosmology: theory - galaxies: formation - galaxies: evolution - quasars: general 



1. INTRODUCTION 

There is a well determined set of cosmological parameters 
h = tfo/lOOkms -1 Mpc -1 = 0.70 ± 0.04, Vt M = 0.27 ± 0.04, 
ft A = 0.73 ±0.05, f = 13.7 ±0.2 Gyr, cr 8 = 0.84 ±0.04 emerg- 
ing from a number of observations (the concordance cos- 
mology, see Spergel et al. 2003). Also the cosmic density 
of baryons ilb = 0.044 ± 0.004 has been very precisely de- 
termined through both the power spectrum of Cosmic Mi- 
crowave Background anisotropies and measurements of the 
primordial abundance of light elements (Cyburt et al. 2001; 
Olive 2002). An important complementary information is 
that the density of baryons residing in virialized structures 
and associated to detectable emissions is much smaller than 
Qfe. In fact, traced-by-light baryons in stars and in cold 
gaseous disks of galaxies and in hot gas in clusters amount 
to a fl bMm w (3-4) x 10" 3 < O.lOfc (Persic & Salucci 1992, 
Fugukita et al. 1998, Fukugita & Peebles 2004). On the other 
hand, in rich galaxy clusters the ratio between the mass of 
the dark matter (DM) component and the mass of the baryon 
component, mainly in the hot intergalactic gas, practically 
matches the cosmic ratio Clju/^b- 

The circumstance that O/, is a factor of about 10 larger than 
^&,ium puts forth both an observational and a theoretical prob- 
lem. On one side, observations are needed to detect and locate 
these "missing" baryons (see for a review Stocke, Penton & 
Shull 2003). On the theoretical side, galaxy formation mod- 
els have to cope with the small amount of baryons presently 
in gas and stars inside galaxies. No doubt that feedback from 
stars and AGNs played a relevant role in unbinding large 
amounts of gas and eventually removing them from the host 
DM halo (see, e.g., Dekel & Silk 1986; Silk & Rees 1998; 
Granato et al. 2001, 2004; Hopkins et al. 2005; Lapi, Cava- 
liere & Menci 2005), but we need to get a detailed quantitative 
understanding of these processes, which are crucial to com- 
prehend galaxy formation. The stellar feedback is expected 
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to depend on the star formation history and the total energy 
released to the gas ultimately depends on the total mass of 
formed stars and on the present-day galaxy luminosity. Cor- 
respondingly, the total energy injected by the AGN feedback 
is ultimately related to the final black hole (BH) mass. The 
fraction of the gas removed by feedbacks is expected to de- 
pend also on the binding energy of the gas itself, which is 
determined by the galaxy virial mass and by its density dis- 
tribution. Therefore the relationships between the galaxy halo 
mass and the galaxy luminosity, the stellar and baryonic mass, 
the mass of the central BH, are expected to give extremely 
useful information on the process of galaxy formation and 
evolution. An additional relevant piece of information is the 
link between the galaxy halo mass and the velocity disper- 
sion of the old stellar component. This paper is devoted to a 
statistical study of these relations. 

The Halo Occupation Distribution (see, e.g., Kauffmann, 
Nusser & Steinmetz 1997; Peacock & Smith 2000; Berlind & 
Weinberg 2002; Magliocchetti & Porciani 2003), which spec- 
ifies the probability P(N\M) that a halo of mass M is hosting N 
galaxies, is a helpful statistics to establish the link between the 
host DM halo mass and the observed galaxy luminosity. An 
additional tool to explore this link is the formalism of the Con- 
ditional Luminosity Function (Yang, Mo & van den Bosch 
2003), which describes how many galaxies of given luminos- 
ity reside in a halo of given mass. Following this approach, 
Yang et al. (2003) investigated the relation between halo mass 
and luminosity. However, particularly for high halo masses, 
Mf, > 10 I3 M Q , both methods give information on galaxy sys- 
tems, more than on large galaxies. 

Marinoni & Hudson (2002) and Vale & Ostriker (2004) 
suggest that a helpful starting point can be the simple hypoth- 
esis that there is a one-to-one, monotonic correspondence be- 
tween halo mass and resident galaxy luminosity. Then, by 
equating the integral number density of galaxies as a func- 
tion of their luminosity and stellar mass to the number den- 
sity of galaxy halos, one gets a statistical estimate of the DM 
halo mass associated to galaxies of fixed luminosity or fixed 
baryon/stellar mass. However, a major problem of the method 
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is the estimate of the mass function of halos hosting one single 
galaxy, the Galaxy Halo Mass Function (GHMF). To solve the 
problem, in this paper we use an empirical approach, which 
takes into account the results of numerical simulations (see, 
e.g., De Lucia et al. 2004 and references therein) on the halo 
occupation distribution by adding to the Halo Mass Function 
(HMF) the contribution of subhalos (Vale & Ostriker 2004; 
van den Bosch, Tormen & Giocoli 2005). At large masses we 
subtract from the HMF the mass function of galaxy groups 
(Girardi & Giuricin 2000; Martinez et al. 2002; Heinamaki et 
al. 2003; Pisani, Ramella & Geller 2003). 

The mass around a galaxy up to a radial distance from its 
center much larger than the characteristic scale of light distri- 
bution can be inferred from detailed X-ray observations (see, 
e.g., O' Sullivan & Ponman 2004). Also the statistical mea- 
surements of the shear induced by weak gravitational lensing 
around galaxies (see Bartelmann, King, & Schneider 2001) 
yield important insights on the halo mass of galaxies (McKay 
et al. 2002, Sheldon et al. 2004). Though these mass esti- 
mates have significant uncertainties and their extrapolations 
to the virial radii are not immediate, they nonetheless provide 
useful reference values to which we compare the outcomes of 
our method. 

The role of stellar and AGN feedback has been discussed 
by several authors (see, e.g., Dekel & Silk 1986; Silk & 
Rees 1998). Granato et al. (2001, 2004) have implemented 
both feedbacks in their model of joint formation of QSOs and 
spheroidal galaxies. More recently the feedbacks have also 
been introduced into hydrodynamical simulations (Springel, 
Di Matteo & Hernquist 2005). One of the purposes of this pa- 
per is to show how the information coming from the relation- 
ships of the galaxy halo mass with measurable galaxy proper- 
ties (such as the stellar, baryonic and central BH masses) can 
shed light on the role and on relative importance of the two 
feedbacks. 

The plan of this work is the following. In §2 we compute 
the galaxy stellar and baryonic mass functions, exploiting the 
luminosity function and the mass-to-light ratio of the stel- 
lar component inferred from kinematical mass modelling of 
galaxies. Then, in § 3, we derive the mass function of galactic 
halos, and, exploiting the relevant galaxy statistics (luminos- 
ity function, galaxy star/baryon mass function, velocity dis- 
persion function) and the galaxy halo mass function, we in- 
vestigate the relationships of the corresponding galaxy prop- 
erties with the halo mass. The relation of the halo mass with 
the mass of the central BH in galactic spheroidal components 
is deduced in §4, by comparing the central black hole mass 
function to the galaxy halo mass function. In § 5 we discuss 
the role of the stellar and AGN feedbacks in shaping the re- 
lationships between stellar and baryonic mass and halo mass. 
§ 6 is devoted to the conclusions. 

2. THE STAR AND BARYON MASS FUNCTIONS OF GALAXIES 

The luminosity function (LF) is a fundamental statistics for 
galaxies. Its present form is the result of physical processes 
involving both baryons and dark matter. In particular, the LFs 
in the range between about 0.1 to several /.im probe the stel- 
lar component. The mass of stars and baryons associated to 
galaxies can be derived coupling the LF with estimates of the 
mass-to-light ratio (MLR) of the stellar and gaseous compo- 
nents. As it is well known, the MLR and the fraction of gas 
depend on galaxy morphology. 

Nakamura et al. (2003) estimated the LF in the r*-band for 
early- and late-type galaxies separately. The separation has 



been done through a light concentration method (see Shankar 
et al. 2004 for a discussion and a comparison with other LF 
estimates). These early- and late-type galaxy LFs are in rea- 
sonable agreement with the LFs of red and blue galaxies, re- 
spectively, derived by Baldry et al. (2004). 

Since the Nakamura LF of late-type galaxies is well defined 
only for luminosities brighter than M r < -18, we extended it 
to lower luminosity using the results of Zucca et al. (1997) 
and Loveday (1998) and translating them from the fey-band 
to the r*-band using a color (bj-r*) w 0.33, as appropriate 
for star forming irregular galaxies (Fukugita et al. 1995; their 
Table 3, panels (j) and (m) with bj = Bj and r* ~ r'). The 
conversion to solar luminosities has been done taking M,. = 
4.62 (Blanton et al. 2001). The resulting LF is well fit, in the 
range 3 x 1O 7 L <L r <3x 1O U L , by 

(t>(L r )&L r _ / 9.05 x IP" 3 4x 1Q- S \ 
Mpc" 3 ~ I^i-i4g0.0076, + x 4 03 ) ' ( ) 

where x = L r /2A x 1O 8 L . 

The MLR of the stellar component can be derived from 
studies of stellar evolution, with uncertainties associated to 
the poor knowledge of details of the IMF (see, e.g., Fukugita 
et 1998; Bell et al. 2003; Fukugita & Peebles 2004; Baldry et 
al. 2004; Panter, Heavens & Jimenez 2004). A more direct ap- 
proach exploits detailed kinematical and photometric studies 
of galaxies to estimate the amount of mass traced by light and 
the mass of the DM component, taking advantage of their dif- 
ferent distribution inside the galaxies. This method has been 
used by Salucci & Persic (1999), who estimated the stellar and 
gaseous mass as a function of the B-band luminosity for late- 
type galaxies to get the baryon mass Mb w 1.33Mei + M star . 
We have approximated their results for the stellar and the gas 
component, respectively as 

log^=-1.6 + 1.21og^, (2) 
for 1O 7 M < M star < 10 12 M Q , and 

Mm Lb 

log -2 = 1.34 + 0.81 log -S. (3) 

Mq Lq 

in the range 3 x 1O 6 M < M HI < 1O U M . Combining these 
relations with the LF of eq. Q, assuming a Gaussian distri- 
bution around the mean relations with a dispersion of about 
20%, we derived the stellar mass function and the baryonic 
mass function of late-type galaxies. To do that we have taken 
M BQ = 5.48 (Binney & Memfield 1998) and (B-r 1 ) = 0.9 (see 
Table 3, panel (m) of Fukugita et al. 1995). 

A similar approach can be followed for early-type galax- 
ies. For about 9000 such galaxies extracted from the Sloan 
Digital Sky Survey (SDSS) Bernardi et al. (2003) estimated 
the MLR (within the characteristic half luminosity ratio r e , 
in solar units) M/L= 3.6(L/L*}° 15 M /L in the r*-band 
(L* = 2 x 1O 1O L ). These authors derived the mass inside r e 
using the relation M{r e ) = ca 2 r e /G, where a is the central ve- 
locity dispersion, and assuming c = 2. However, the value of 
c depends on the light profile (c ps 2.35 for a de Vaucouleurs 
profile, see Prugniel and Simien 1997) and on the DM distri- 
bution (Borriello et al. 2003). Using c = 2.35 and rescaling the 
zero point of the MLR by Bernardi et al. (2003), we obtain: 
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FIG. 1 . — The Galaxy Stellar Mass Function. The thick solid line shows 
the numerical results, which are almost perfectly matched by the analytic 
fitting formula [eq. j5J]- The dotted line gives the contribution of early-type 
galaxies. The shaded area represents the uncertainty due to the ~ 30% error 
in the mean stellar mass to light ratio. The vertical dot-dashed line shows 
the stellar mass corresponding to M r = —18. The dashed line shows the result 
by Bell et al. (2003). The inset illustrates the average fraction of gas as a 
function of the stellar mass. 



As discussed by Bell et al. (2003) and by Baldry et al. (2004), 
the uncertainties related to the IMF and to the star formation 
history imply an uncertainty of about 30% in the mean value 
of the MLR; the dispersion around the mean relation, eq. 
is of about ~ 20%. By convolving the r*-band luminosity 
function of early-type galaxies of Nakamura et al. (2003) 
with the distribution of MLRs, assuming a Gaussian distribu- 
tion around the mean relation of eq. ®, with a dispersion of 
20%, we estimate the Galaxy Stellar Mass Function (GSMF) 
of E/SQ galaxies, which essentially coincides with the Galaxy 
Baryonic Mass Function (GBMF), since in early-type galax- 
ies the gas gives a negligible contribution to the baryon mass. 

The total GSMF, holding over the mass range 10 8 M Q < 
Mstar ^ 10 12 M Q , is shown by the solid line in Fig. [J where the 
shaded area corresponds to the 30% uncertainty in the mean 
MLR. This is a safe range to determine the GSMF as at lower 
masses the uncertainty in the LF grows while the increase 
in the total stellar mass density including M star < 10 8 M Q is 
rather small, < 5%. The upturn at M stal <3x 10 8 M Q cor- 
responds to the appearance of the dwarf galaxy population, 
represented by the second term at r.h.s. of eq. Q. How- 
ever the contribution in stellar mass density in the range 
10 8 M Q < M star < 10 9 M Q is just - 3% of the total. In the 
inset of Fig.[^we have displayed the gas fraction as a func- 
tion of the stellar mass. A very accurate analytical represen- 
tation (actually indistinguishable from the solid line showing 
the numerical results) is provided by a Schechter function plus 
a power law term: 
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GSMF(M star )dM star 
Mpc" 3 

where x = M star /6 x 10 I0 M Q . The GBMF is easily computed 
by adding the appropriate gas contribution. Recent estimates 
of the GSMF and of the GBMF have been produced by Cole 
et al. (2001), Bell et al. (2003) and Baldry et al. (2004), 
exploiting 2dF, SDSS and Two Micron All Sky Survey data. 
The estimate by Bell et al. (2003) is shown by the dashed 
line in Fig.^ The difference with our estimate is mainly due 
to the difference in adopted MLR. Bell et al. (2003) have 
derived their MLR by fitting the broad band SED with stellar 
population models. The mean MLR adopted here, based on 



kinematic determinations, is a factor of about 1.3 higher at 
high luminosities, while at very low luminosities (L r < 5 x 
1O 8 L0) it is about a factor of 2 lower. However, the flatness 
of the GSMF at small masses conceals the difference in the 
MLR, while at large masses the almost exponential decline of 
the LF amplifies it. It is worth noticing that the determination 
of the MLR of low luminosity objects is hampered by many 
effects related to the episodic star formation history, to the 
presence of dust, to the irregularity of their shapes, and to the 
DM predominance. 

Cole et al. (2001) presented estimates of the GSMF for two 
choices of the Initial Mass Function (IMF): Salpeter's (1955) 
and Kennicutt's (1983). For a Salpeter IMF their GSMF is 
consistent with ours, within our estimated uncertainties: the 
main difference is a ~ 30% excess for log(M stal -/M0) < 10.5. 
For a Kennicutt IMF, their GSMF drops for M stal about 0.2 
dex lower than that by Bell et al. (2003). The estimate by 
Baldry et al. (2004) is close to that by Bell et al. (2003), as 
expected since they exploit very similar LFs and MLRs. All in 
all, methods based on kinematic measurements and on stellar 
population synthesis yield GSMFs and GBMFs in reasonable 
agreement, and establish a sound confidence interval. 

Integrating the GSMF for M star > 10 8 M Q , the mass density 
parameter of baryons condensed in stars associated to late- 
type galaxies is found to be 



^ tar (kin) = (1.8±0.4)x 10" 



(6) 

where the label "kin" indicates that the stellar mass of galaxies 
has been estimated using kinematic data. The corresponding 
neutral gas density amounts to ~ 20% of this value and it is 
concentrated in late-type, low mass systems with M star < 5 x 
10 9 Mq. Here and in the following eqs. Q and (|8j, the errors 
reflect the uncertainties on the GSMF. 

The star density parameter associated to early-type galaxies 
amounts to 



r»f tar (kin): 



(1.8 ±0.6) x 10" 



(7) 



It is well known that in early-types the amount of cold gas is 
negligible. Therefore, the overall local stellar mass density 
in galaxies with stellar masses in the range 10 8 M Q < Mt, < 
10 12 M M is 



^ ar (kin) = (3.6±0.7)x 10 -3 . 



(8) 



This value is in good agreement with the recent estimates ob- 
tained through spectro-photometric galaxy models (Bell et al. 
2003; Fukugita 2004; Fukugita & Peebles 2004). The contri- 
bution of the cold gas to the baryon density in galaxies is only 
~ 8%, and thus flf « 1.080^. This result confirms the well 
known conclusion that only a small fraction, < 10%, of the 
cosmic baryons is today in stars and cold gas within galax- 
ies. It is worth noticing that the star formation rate integrated 
over the cosmic history matches the overall local mass den- 
sity in stars (see Nagamine et al. 2004). This mass density 
has been accumulated at high redshifts, z > 1, for early-type 
galaxies and later on for late-types, as indicated by their re- 
spective stellar populations. 

By subtracting from the cosmic matter density the contri- 
bution of baryons (Q b < 0.044) and that of dark matter (DM) 
in groups and clusters of galaxies (£1% M ~ 0.012; Reiprich 
& Bohringer 2002), we obtain the DM mass density associ- 
ated with galaxies £l£ M w 0.21, in excellent agreement with 
the determination by Fukugita & Peebles (2004). The aver- 
age DM-to-baryon (essentially stars) mass fraction in galaxies 
turns out to be around 60. This value must be compared with 
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the cosmological ratio R cosm = flnnf/tib ~ 6. In fact, in rich 
galaxy clusters the baryon mass, mostly in form of diffuse gas, 
and the DM halo mass have a ratio consistent with the "cos- 
mic" DM to baryon ratio (see, e.g., Ettori, Tozzi & Rosati 
2003). This evidences that the baryon fraction in galaxies de- 
creases on average by a factor of about 10 relative to the initial 
value, due to a number of astrophysical processes associated 
to the formation of these objects. In the following we will use 
the cosmic fraction / cosm = 1 /R cosm ~ 0.17. 

3. THE GALAXY HALO MASS FUNCTION AND THE L, M ST ar AND 
a VS. HALO MASS RELATIONS 

In order to investigate the relationships between the stellar 
and baryonic mass and the halo mass in galaxies in a one- 
to-one correspondence, the statistics of halos containing one 
single galaxy, the Galaxy Halo Mass Function (GHMF), has 
to be estimated. The overall HMF as found by numerical sim- 
ulations (see, e.g., Jenkins et al. 2001; Springel et al. 2005) 
is well described by the Press & Schechter (1974) formula as 
modified by Sheth & Tormen (1999). However, in order to 
compute the GHMF, we have to deal with the problem of the 
halo occupation distribution (HOD; Peacock & Smith 2000; 
Berlind et al. 2003; Magliocchetti & Porciani 2003; Kravtsov 
et al. 2004; Abazajian et al. 2005). Two effects need to 
be taken into account: (i) the number of galaxies is actually 
larger than the number of DM halos, since a halo may con- 
tain a number of sub-halos, each hosting a galaxy, and (ii) the 
probability that very massive halos (M h > 10 12 -10 13 M Q ) host 
a single giant galaxy drops rapidly with increasing halo mass. 

To account for the effect (i) we use the results by Vale & 
Ostriker [2004; see their eqs. (1) and (3)] and we add to the 
HMF the subhalo mass function they have derived; we have 
checked that this procedure does not alter substantially the 
overall mass density in the galactic range. The subhalo MF 
by van den Bosch et al. (2005) is extremely close to the Vale 
& Ostriker one over the mass range of interest here, and gives 
essentially indistinguishable results. To account for effect (ii) 
we subtract from the HMF the halo mass function of galaxy 
groups and clusters. Estimates of the latter obtained by dif- 
ferent groups (Girardi & Giuricin 2000; Martinez et al. 2002; 
Heinamaki et al. 2003; Pisani et al. 2003; Eke et al. 2006) 
are in reasonable agreement for Mi, >5x 10 12 M Q . At lower 
masses the galaxy group mass function is very uncertain, and 
the recent study by Eke et al. (2006) finds a larger abun- 
dance of low luminosity groups than previously reported from 
smaller samples. On the other hand, from Fig. 8 (plus Fig. 15) 
of Eke et al. (2006), it looks plausible that galaxies dominate 
the MF for M h < 5 x 1O 12 M . We stress, however, that, as 
discussed in the following, in the mass range considered here 
our results are only weakly sensitive to whether or not the 
galaxy group mass function is subtracted from the HMF. 

The GHMF obtained subtracting from the HMF the group 
and cluster mass function estimated by Martinez et al. (2002) 
is shown in Fig. |2] and is well fit, in the range 11 < 
logM/,/M Q < 13.2, by a Schechter function 



(M h 



GHMF(M A ) dM,, = -r 

M V M 



-m/M) 



dM h 



(9) 



witha = -1.84,M=1.12x 10 13 M© and6> = 3.1 x lO^Mpc" 3 . 
The fall off at high masses (where early-type galaxies domi- 
nate) mirrors the increasing probability of multiple occupa- 
tion of mass halos found by Magliocchetti & Porciani (2003) 
forM > 3 x 10 13 M o (see also Zehavi et al. 2005). Weak lens- 



O 



Galac:ic Halo Mass Function 
Halo + SubHalo Mass Function 
Martinez et al, (2002) 
Group Halo Mass Function (Fit) 



10.5 11.0 11.5 12.0 12.5 13.0 13.5 14.0 
Log M„ [M e ] 

FIG. 2. — Galactic halo mass function. The heavy solid line shows the 
numerical results, obtained as described in Sect.0 which are very accurately 
reproduced by the analytical approximation [eq. J9j]. We have plotted the 
results by Martinez et al. (2002) only for the halo mass range of interest here. 
Beyond M/, fa 2 X 1O 13 M0, the halo mass function is fully accounted for by 
groups and clusters of galaxies; as a consequence the GHMF goes to zero 
(shown as a dotted line beyond this limit). 



ing measurements also suggest an upper galaxy mass limit 
M max < 3 x 10 13 M Q (Kochanek & White 2001). 

If two galaxy properties, q and p, obey a monotonic rela- 
tionship we can write 



Hp)^rdq=^(q)dq, 
dq 



(10) 



where ^(q) is the number density of galaxies with measured 
property between q and q + dq and <£>(/?) is the corresponding 
number density for the variable p. The solution is based on 
a numerical scheme that imposes that the number of galaxies 
with q above a certain value q must be equal to the number of 
galaxy halos with p above p (see, e.g., Marinoni & Hudson 
2002; Vale & Ostriker 2004), i.e., 



<S>(p)dp -- 



Wq)dq 



(11) 



In the following p = M h and $(p) = GHMF(M/,), while the 
variable q will be, in turn, the luminosity, the stellar mass, 
the velocity dispersion, and the central BH mass. It is worth 
noticing that in this way we establish a direct link between 
the specific galaxy property and the halo mass without any 
assumption or extrapolation of the DM density profile. 

The monotonicity assumption is obviously critical for the 
applicability of the present approach. However direct ev- 
idence of monotonic relationships between several pairs of 
these quantities has been reported (Haring & Rix 2004; Fer- 
rarese 2002; Baes et al. 2003; Tremaine et al. 2002), and 
additional data supporting the relationships derived here are 
discussed in the following. 

We also implicitly assume that all galactic-size halos con- 
tain a visible galaxy. This assumption underlies all ma- 
jor semi-analytic models for galaxy formation, including the 
one by Granato et al. (2004), that we adopt as our refer- 
ence model. The successes of this model in reproducing 
the redshift-dependent galaxy luminosity function in different 
wavebands provide strong support to this view. 

The result for the stellar mass, obtained setting ^f(q) = 
GSMF(M star ) [eq. @] is plotted in Fig. We find that its 
relationship with the host halo mass is well approximated by 
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FIG. 3. — (a) Mass in stars versus halo mass. The thick solid line shows 
the numerical results while the dot-dashed line (difficult t o di stinguish from 
the solid line) represents the analytic fitting formula [eq. 1121 1. The dashed 
line has been obtained using the GSMF by Bell et al. (2003); the barred area 
represents the uncertainty associated to the mass to light ratio; the shaded 
area illustrates the result at z = 1.75, based on data by Fontana et al. (2004). 
(b) r*-band luminosity as a function of halo mass. Again, the thic k so lid line 
shows the numerical results and the analytic fitting formula [eq. 1131 1 is in- 
distinguishable from it. The barred area represents the uncertainty associated 
to the LF (dominant at the low-mass end) and to the GHMF (dominant at the 
high-mass end); the thin dashed line is the numerical result obtained without 
removing the galaxy groups and clusters from the HMF; the dot-dashed line 
is the result by Vale & Ostriker (2004); the shaded region shows the result by 
Kleinheinrich et al. (2004); the arrow is from O'Sullivan & Ponman (2004); 
the diamond is from Hoekstra et al. (2004); the star is from Guzik & Seljak 
(2002). (c) cr-M;, relation. As before, the thick solid line shows the numer- 
ical results while the dot-dashed line shows the analytic formula [eq. 1151 1 
holding in the large M;, limit. The barred area represents the uncertainty. 



The calculations for the baryonic mass is strictly analogous. 
Setting ^{q) = <fi(L r ) [eq. Q] we also derived the approxi- 
mated behavior of the luminosity as a function of the halo 

mass 



L r , 



1.2x10'%, W3xlO»tf ) 



2.65 



l+(M h /3 x lO 11 *^) 2 - 00 
and of the halo mass as a function of luminosity 



Mi, 



3 x 10 n M ( 



© 



1.3 x 10 10 L 







0.35 



(13) 



1 .65 



1.3 x 10 10 L 







(14) 

Both stellar mass and luminosity exhibit a double power law 
dependence on halo mass with a break around M/, b rea k ~3x 
10 n M Q , corresponding to a luminosity L r ,break ~6x 10 9 L Q . 

The derivation of the L r -M/, relation is quite sensitive to un- 
certainties in the LF and in the GHMF. In the range 10 9 L Q < 
L r ;$ 10 11 Lq the LF is rather precisely known (see the com- 
parison of different LF estimates in Fig. 1 of Shankar et al. 
2004), while the effect of uncertainties in the GHMF becomes 
significant at large masses. On the whole, the L r -Mi, rela- 
tion is quite accurately determined in the range 1O U M0 < 
Mi, < 10 13 M Q (Fig.[3J>). At low luminosities the errors on 
the LF rapidly increase, becoming a factor of about 2 for 
Lr S 3 x 10 8 L©. In order to illustrate the consequences on 
the L r -Mi, relation, we have computed it using the 1 a upper 
and lower boundaries of the LF by Nakamura et al. (2003). In 
the former case, the low-luminosity portion of the L r -Mh rela- 
tion flattens from a slope of ~ 2.6 to ~ 1 .9; in the latter case, it 
decreases almost exponentially. Therefore, the extrapolation 
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4. — Halo to stellar mass ratio as a function of the stellar mass. The 



solid line is the result of numerical calculations using e 

qJTT} with <S>(p) : 

GHMF(M^) and *(g) = GSMF(M star ), as given by [eq. gj]; the barred area 
represents the uncertainty associated to the mass to light ratio. The dashed 
line has been obtained using the GSMF by Bell et al. (2003). The Cole et 
al. (2001) GSMF with a Salpeter (1955) IMF yields results very close to 
ours. The dark shaded area represents the data on giant elliptical galaxies by 
Gerhard et al. (2001); the light shaded area represents data on spiral galaxies 
by Persic, Salucci & Stel (1996) and Salucci & Burkert (2000). 



of the above relationships to L r <C 10 9 L© and, correspond- 
ingly, to Mi, <C 10 n M Q must be taken cautiously. Of course, 
the same conclusion holds for relationships involving other 
statistics of galaxies and of their host halos. This emphasizes 
the need of a precise determination of the low luminosity end 
oftheLF. 

In Fig. [5J> our estimate of the L r vs. Mi, relation is com- 
pared with observational results for galaxies whose M/, could 
be derived based on two different methods: (i) an X-ray- 
based mass model of the isolated elliptical galaxy NGC 4555 
(O'Sullivan & Ponman 2004), in which the gravitational po- 
tential is known up to about 1/8 of the virial radius (the mass 
within this radius is shown as a lower limit in Fig. |3j>); (ii) 
weak-lensing observations that provide the shear field around 
a number of galaxies of average luminosity L, from which 
it is possible to infer the projected mass density and even- 
tually to extrapolate the virial mass by assuming a DM pro- 
file (Guzik & Seljak 2002; Kleinheinrich et al. 2004; Hoek- 
stra et al. 2004). In Fig. [3J> we also show for compari- 
son the L r vs. Mi, relation obtained from eq. (II Oi setting 
<l>(p) = HMF(Mi,). It is worth noticing that if the group and 
cluster MF is not subtracted from the HMF, our results vary 
only at high masses, and the changes do not exceed 0.2 dex 
uptoM/, ~2 x 10 13 M Q . 

As a further check, our estimate of the ratio M/,/M sta , as a 
function of M star , is compared in Fig.|4]with estimates derived 
by extending to the virial radius the inner mass models of a 
number of giant ellipticals (Gerhard et al. 2001) and spirals 
(Persic, Salucci & Stel 1996; Salucci & Burkert 2000). We 
stress that these results require extrapolations to the virial ra- 
dius of the density profile, assumed to have the Navarro, Frenk 
& White (1996) shape, while our estimate does not need any 
assumption on DM density profile. It is apparent that these 
independent results are in nice agreement with our findings. 

The dependence of the luminosity on the halo mass has 
been investigated also by Vale and Ostriker (2004); their re- 
sult is also shown in Fig.[3J>. They exploited the 2dF galaxy 
luminosity function in the /by-band (Norberg et al. 2002), ex- 
trapolating it beyond the range of magnitudes where it was 
defined. The difference in the L r -Mf, relation between our es- 
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timate and theirs is due to the steeper slope of the LF adopted 
by them. There is also a small difference in the normalization 
of the LF, but this is of minor importance. 

At high masses the direct comparison of the galaxy LF to 
the halo plus subhalo number density [cf. their eq. (9)] results 
in a slight flattening of the relation (dot-dashed curve for Vale 
& Ostriker, dashed curve for our calculations). As Vale and 
Ostriker (2004) pointed out, in this way the mass term refers 
to the entire halo hosting the group or the cluster and not to 
just the galaxy halo. 

Guzik & Seljak (2002) modelled the galaxy-galaxy lensing 
trying to separate the central galactic contributions from con- 
tributions of the surrounding groups and clusters. Their model 
applied to the SDSS data on galaxy lensing yields Mi,/L r w 
50Mq/Lq at the characteristic luminosity 3 x 10 10 L© for 
early-type galaxies, in keeping with our results. They also 
found a luminosity dependence M;, oc L 1 4±0 2 , compatible 
with the high-luminosity slope of eq. J14t . However, we find 
that, at low luminosities, the slope significantly flattens to- 
ward a dependence M/, oc L 035 , while Guzik & Seljak (2002) 
assume a single power law relation. 

Van den Bosch et al. (2003, 2005) computed the con- 
ditional luminosity function of early- and late-type galax- 
ies, a statistics linking the distribution of galaxies to that 
of the DM. They concluded that the MLR has a minimum 
M h /L r ~45-7OM Q /L at M/, ~ 2 — 4 x 10 h Mq, consistent 
with eq. ( fT4l . 

Marinoni & Hudson (2002) investigated the MLRs of the 
virialized systems, which include galaxies, groups and clus- 
ters. By comparing their luminosity function to the ACDM 
halo mass function, they concluded that the MLR has a broad 
minimum around Lj«3x 10 10 Lq. The slopes at low and 
high luminosity are -0.5 and +0.5, respectively. Our slope 
is similar at low masses, where practically all virialized sys- 
tems are galaxies and thus the comparison is meaningful. The 
studies by Eke et al. (2004, 2006) of the variation of the MLR 
with size of galaxy groups is fully consistent with a minimum 
at approximately the same luminosity. 

By comparing the HMF and the LF, as we have done for lo- 
cal galaxies, it is possible to infer the M stal — M/, relation even 
at substantial redshifts. For the GSMF we use a simple linear 
fit [log(^(M star )/Mpc 3 ) = -1.71og(M star /M Q )+16.1, holding 
for 11 < log(M stal /M Q ) < 12] to the data by Fontana et al. 
(2004) at z = 1.75, and we approximate the GHMF with the 
HMF computed at the same redshift. The result, shown by the 
shaded area in Fig.|3j?, has to be taken as an upper limit since 
we have neglected the contribution of galaxy groups to the 
HMF. Clearly, our estimate becomes increasingly uncertain 
as we approach the upper limit of the interval where GSMF 
is observationally estimated; this is reflected in the increased 
width of the shaded area. Nevertheless, the M stal — M/, relation 
at z = 1 .75 turns out to be quite close to the local one, indicat- 
ing that for large galaxies the M stal — M/, relation was already in 
place at redshift z > 1 in line with the anti-hierarchical baryon 
collapse scenario developed by Granato et al. (2001; 2004). 

Sheth et al. (2003) estimated the Velocity Dispersion Func- 
tion (VDF) of spheroidal galaxies using a sample drawn from 
the SDSS survey and have built a simple model for the con- 
tribution to the VDF of the bulges of spirals, which dominate 
at low velocity dispersions. Their estimate covers the range 
80 km s" 1 < a < 400 km s" 1 . Comparing the global VDF (in- 
cluding both early and late type galaxies, as shown in Fig. 6 
of Sheth et al.), and the GHMF with the same technique pre- 
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FIG. 5. — Fraction of primordial gas turned into stars as a fu nctio n of 
halo mass. The solid line has been obtained numerically from eq. Jll| with 
*(p) = GHMF(M A ) and f(q) = GSMF(M st!11 ■), as given by [eq. gjj; the 
barred area represents the uncertainty associated to the mass to light ratio. 
The dashed line has been obtained using the GSMF by Bell et al. (2003). We 
have set / C osm = 1/6. 



sented above [cf. eqs. JlOi and ( II 1H . we can derive the cr-Mh 
relationship shown in panel (c) of Fig. [5] ForM/, >6.3 x 10 11 , 
the relationship is accurately represented by a simple power- 
law: 

1 /3 

crwllOkms" 1 ! — h — ) , (15) 

V6.3xl0 n M Q y ' V ; 

while it steepens for lower halo masses. The uncertain- 
ties strongly increase for a < 80 km s" 1 , corresponding to 
Mh < lO n M . Note that these relationships must break down 
in the low-er (and low-M/,) regime. In fact, the close match 
found by Cirasuolo et al. (2005) between the VDF and the 
virial velocity function (derived from the halo mass distri- 
bution function integrated over redshift) indicates that halos 
more massive than M/, ~ 10 11 M Q are generally associated to 
spheroidal galaxies or to later type galaxies with bulge veloc- 
ity dispersions a > 80 km s" 1 . On the other hand, the fraction 
of galactic halos associated to essentially bulgeless late-type 
galaxies must increase with decreasing M/,, so that the integral 
of the VDF deviates from that of the GHMF and eq. il It no 
longer applies. 

Figure [5] shows the ratio of the mass in stars to the ini- 
tial baryon mass associated with each halo, assumed to be 
Mt.i = fcosmM/,, as a function of M/,. It illustrates the "inef- 
ficiency" of galaxies, especially of those of low halo mass, 
in retaining baryons. As discussed in §|5J the shape of the 
Msxzt/ fcosmMh can be understood in terms of feedbacks: at 
low masses the SN feedback is very efficient in removing the 
gas, thus quenching the star formation; moving toward higher 
masses (for log(MfiH) > 7.5, corresponding to log(M/, > 12)) 
the AGN feedback becomes more and more powerful, to the 
point of sweeping out most of the initial baryons. 

Our result is at odds with the claim by Guzik & Seljak 
(2002) of a high efficiency, up to 75%, in turning primor- 
dial gas into stars. However, the claim is based on a MLR 
Mh/Li w HIiMq/Lq in the i' band for late-type galaxies, a 
factor of 3 lower than the value found for the early-type ones. 
As the authors themselves point out, the statistical signifi- 
cance of this result is marginal, due to the weak lensing signal 
for the fainter late-type galaxy sample. The GSMFs by Cole 
et al. (2001) and Bell et al. (2003) yield similar efficiencies, 
which are very close to our estimates for relatively low halo 
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masses, but lower for large masses, yet within the estimated 
uncertainties. 

4. BLACK HOLE VS. HALO MASS 

The relation between the central supermassive black hole 
and the halo mass M/, is relevant in the framework of theories 
for the origin and evolution of both galaxies and AGNs (Silk 
& Rees 1998; Monaco, Salucci & Danese 2000; Granato et al. 
2001; Ferrarese 2002; Granato et al. 2004). To constrain such 
relation we adopted the procedure presented in the previous 
Section [eqs. dlOl and il 1H . replacing the function ty(q) with 
the local BH mass function (Shankar et al. 2004). We assume 
that each galactic halo hosts just one supermassive BH. Our 
result is shown in Fig.|6] where the barred area illustrates the 
errors due to the observational uncertainties on the BH mass 
function, as estimated by Shankar et al. (2004). We find good 
agreement, within the estimated uncertainties, with the pre- 
dictions of the Granato et al. (2004) model. 

The relationship can be approximated by 



3.95 



Mr 



: 6 x 10 b M ( 







(M h /2.2 x 10 n M Q ) 
\ + {M h /2.2 x 1O U M ) 2 



(16) 



Again a double power law with a break at Mi, ~3x 10 11 M Q 
is a very good representation of our results. At the high mass 
end, the BH mass is nearly proportional to the halo mass 
(Mrh oc M^ 25 ), while at low masses the relation steepens sub- 
stantially (M BH oc Ml 95 ). 

In Fig. [6] we also compare our estimate of the Mbh-M/, 
relation with that of Ferrarese (2002), who first investigated 
this issue from an observational point of view. She derived a 
power-law relationship between the bulge velocity dispersion 
and the maximum circular velocity, v c , for a sample of spiral 
and elliptical galaxies spanning the range 100 < v c < 300 km 
s , and combined it with the relationship between v c and the 
virial velocity, Vvir based on the numerical simulations by Bul- 
lock et al. (2001) and with the M/,-v v i r relationship given by 
the ACDM model of the latter authors for a virialization red- 
shift Zvir ~ 0, to obtain a M/,-cr relation. Coupling it with one 
version of the observed BH mass vs. stellar velocity disper- 
sion relationship (Mrh oc cr 4,58 ) she obtained Mbh ocM ; ", with 
a = 1.65-1.82. Baes et al. (2003) with the same method, but 
assuming Mbh oc cr 4 02 and with new velocity dispersion mea- 
surements of spiral galaxies with extended rotation curves, 
yielding a slightly different v c -a relation, found M BH ocM ; ' 11 . 
It is apparent from Fig.[6]that our result differs substantially in 
normalization, while the high-mass slope is remarkably close 
to that obtained by Baes et al. (2003). It should be noted that 
the M/,-V v i r relation depends on the virialization redshift. For 
z v i r ~ 3 [the median virialization redshift for galaxies with ve- 
locity dispersions in the range probed by Ferrarese and Baes 
et al., according to the analysis by Cirasuolo et al. (2005; see 
their Fig. 1)], its coefficient would be a factor of ~ 4.25 lower 
than that used by Ferrarese (2002) and Baes et al. (2003) and 
the coefficients of the Mbh oc M/, relations would be larger by 
a factor of ~ 5.6 in the case of eq. (6) of Ferrarese(2002) or 
of ~ 4 in the case of Baes et al. (2003), bringing them much 
closer to ours. In fact, as suggested by Loeb & Peebles (2003) 
and shown in details by Cirasuolo et al. (2005), the velocity 
dispersion of the old stellar population (whose mass is related 
to the central BH mass) is closely linked to halo mass and 
characteristic velocity at the virialization redshift. 

As a consistency test, we have combined the Mbh-M/, re- 
lation, shown in the upper panel of Fig. [6] with the cr-M/, 
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FIG. 6. — Upper panel: super-massive BH m ass v ersus halo mass. Thesolid 
line has been obtained numerically from eq. 1111 with $(/>) = GHMF(M/,) 
and &(q) is the local super-massive BH mass function estimated by Shankar 
et al. (2004). The dashed line is the relation by Ferrarese et al. [2002; 
her eq. (6)]; the dot-dashed line is the relation by Baes et al. (2003). The 
shaded area represents the prediction of the Granato et al. (2004) model. 
Lower panel: Mbh-"" relation obtained combining the M^u-Mf, relation (up- 
per panel) with the cr-M;, relation (panel (c) of Fig. |5J. The data are from 
Ferrarese & Ford (2005). In both panels the barred area reflects the uncer- 
tainty associated to the BH Mass Function. 



relation, shown in panel (c) of Fig. [3] to obtain the Mbh-c re- 
lation (see Ferrarese_& Ford 2005 for a review), which turns 
out to be consistent with the data, as shown by the lower panel 
in Fig.|6] 

5. FEEDBACK FROM STARS AND AGNS 

The dependence of the star and BH masses on the halo 
mass found in the previous Section, suggests that different 
physical mechanisms are controlling the star formation and 
the BH growth above and below M/,.b re ak ~3 x 10 u Mq, cor- 
responding to M star ~ 1.2x 10 10 M Q after ecL (fl2l . and to 

L r ~5x 1O 9 L (or M, 19.6) after eq. O. It is worth 

noticing that the analysis of a huge sample of galaxies drawn 
from the SDSS shows that around M star w 2-3 x 10 10 M o and 
M r ~ -19.8 there is a sort of transition in the structure and 
stellar ages of galaxies (Kauffmann et al. 2003; Baldry et al. 
2004). 

The efficiency of star formation within galactic halos of dif- 
ferent masses is the result of several processes. The most im- 
portant are: (i) the cooling of the primordial gas within the 
virialized halos (White & Rees 1978); (ii) the injection of 
large amounts of energy into the ISM by supernova explo- 
sions (Dekel & Silk 1986; White & Frenk 1991; Granato et 
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FIG. 7. — Specific energy feedback from stars and AGNs compared to 
the baryon specific binding energy within the host halo, as a function of the 
halo mass. The solid line refers to our estimate of the SN feedback while the 
dashed line is the SN feedback obtained using the GSMF of Bell et al. (2003). 
The dot-dashed line is our estimate of the AGN feedback. The shaded area 
shows the specific binding energy of the gas in the DM potential well, for 
virialization redshifts < z v ; r < 5. 



al. 2001; Romano et al. 2002) and by the central quasar (Silk 
& Rees 1998; Granato et al. 2001, 2004; Lapi et al. 2005). 
All these processes have been implemented in the model of 
Granato et al. (2004). A simplified, analytical rendition of 
this model is presented in the Appendix A. 

The impact of stellar and AGN feedback is illustrated by 
Fig. The binding energy of baryons in the DM potential 



well per unit baryonic mass as a function of the halo mass [cf. 
eq. (IA4H for < z v ir < 5 is shown by the shaded area. To 
compute the overall energy per unit baryonic mass injected in 
the gas by supernovae { E stM ) a nd by AGNs (E A gn), we have 
exploited eqs. ( IA9> and (IA14> . respectively, where M star and 
Mbh as functions of the halo mass are given by eqs. (I12> and 
d!6i . respectively. Then we divided the overall energies by 
the initial baryon mass Afj j = / cosm M;,. Figure shows that 
for large masses the gas can be efficiently removed by the 
AGN feedback which overwhelms the binding energy. For 
small masses the supernova feedback dominates but appears 
to be insufficient to remove the gas associated to the host halo, 
due to the above mentioned problem that the observed M star - 
Mf, relation inferred from the data exhibits a too steep low- 
mass slope. The flattest slope allowed by the data, discussed 
in §|3] would largely alleviate, but not completely overcome, 
this problem. 

The relative importance of the two feedbacks depends on 
their efficiency in transferring the available energy to the gas. 
It is interesting that with the efficiencies used in Fig. Qthe 
crossing point is quite close to M/,,b re ak ~ 3 x 10 n M Q . As 
discussed by Granato etal. (2004) and Cirasuolo et al. (2005), 
a more accurate evaluation of the efficiencies can be obtained 
by fitting statistics of galaxies and of AGNs, such as LFs at 
low and high redshift, the Faber and Jackson relation, and the 
local BH mass function. 

A more quantitative insight into the role of the key ingre- 
dients of the model is provided by the analytic calculations 
presented in th e Ap pendix A. So long as the star formation 
rate obeys eq. (IA2> . the mass in stars at a th e present time t, 
assumed to be >> t c , is given, after eq. ( IA8I . by: 



where f surv is the fraction of stars that survive up to now. 

For large halo masses, where the stellar feedback is less ef- 
ficient (a < 1), the quantity 1 —R + a is a slowly decreasing 
function of the halo mass, so that M stal is approximately pro- 
portional to M/,. However, in this case, the fraction of gas 
turned into stars is controlled by the AGN feedback, which, 
as shown by the full treatment by Granato et al. (2004), for 
M/, >3x 10 n M© expels an approximately constant fraction 
of the initial gas, thus preserving the approximate proportion- 
ality between M star and Mi,, in agreement with eq. dl 2i . 

The effective optical depth, which rules the flow o f the cold 
gas into the reservoir around the BH [cf. eq. (A101], is large 
(r > 1) for large galaxies, implying, after eq. (IA1 It . Mbh ~ 
1.2 x 10~ 3 M star . As a consequence, in the high mass limit, 
the BH mass must have a dependence on the halo mass very 
similar to that of the stellar mass, in agreement with eq. d!6t . 

For Mi, <C 10 12 M Q the dominant term in the denominator 
in the right-hand side of eq. (I17> is the effective efficiency of 

the SN feedback, a oc M h 2/3 [cf. eq. iAgt l. Therefore we get 



M,.„ x U, Ml ' ■ (18) 
This limiting dependence has been derived theoretically also 
by Dekel & Woo (2003) with similar assumptions. On the 
other hand, such a relation is significantly flatter than that in- 
ferred from the data [cf. eq. dl2l and Fig.Q. Possible expla- 
nations may be that in less massive halos the initial baryon 
fraction is lower by effect of reheating of the intergalactic 
medium, hindering the infall of baryons into the shallower po- 
tential wells, or that the SN efficiency in removing the gas is 
higher. However, the difference must not be overemphasized, 
in view of the large uncertainties on the shape of the M star -M;, 
and M/,-L r relations at low masses/luminosities, induced by 
our poor knowledge of the low-luminosity portion of the LF. 
As discussed in §13 the data are consistent also with a flatter 
relation (M star oc M ; ' 9 ). 

Since in the mass rang e My, < 10 u Mq the optical depth is 
small (r <C 1), from eqs. JA1 1> and (IA12i we obtain: 

M BH n(M star r«Mf . (19) 

Thus this simple model predicts that the low mass slope of 
the Mbu-Mi, relation is steeper than that of the M star -M/, re- 
lation, because of the decrease of the optical depth with mass 

2/3 

r oc M h , entailing a lower capability of feeding the reservoir 
around the BH. Interestingly, a steepening by approximately 
this amount is also found from our analysis of observational 
data [cf. eqs. (O and (tl6tl. 

6. CONCLUSIONS 

We computed the stellar and baryon mass function in galax- 
ies exploiting M/L ratios for stars and gas derived from 
galaxy kinematics. The results turn out to be in agreement 
with previous analyses based on stellar population models. 
The total baryonic mass density in galactic structures amounts 
to ft£ w (3.6 ± 0.7) x 10" 3 , of which - 40% resides in late- 
type galaxies. This result confirms the well-known conclusion 
that only a fraction < 10% of the cosmic baryons are presently 
in stars and cold gas within galaxies. 

The present-day galaxy halo mass function, i.e., the num- 
ber density of halos of mass M/, containing a single bary- 
onic core, has been estimated by adding the subhalos to the 
halo mass function and by subtracting the contributions from 
galaxy groups and clusters. Such subtraction, which is re- 
quired to single out galactic halos, has however a minor effect 
over the mass range of interest here. 



M star OC / Sl 



l-R + a 



(17) 
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Approximated analytic relationships between the halo 
mass, Mh, the mass in stars, M stal -, and the r*-band lumi- 
nosity, L r , have been obtained from the functional equations 
<fr[> M h (q)] = *(> q), where $(> M h ) is the number den- 
sity of galactic halos larger that Mh and >J/(> q) is the num- 
ber density of galaxies with either stellar mass greater than 
M star or luminosity greater than L r . The results are in good 
agreement with M/,/L r ratios inferred through X-ray map- 
ping of the gravitational potential and through gravitational 
lensing. Both relations exhibit a double-power law shape 
with a break around M/,b rea k ~ 3 x 1O U M0, corresponding 
to M star ,break ~ 1.2 x 10 10 M Q and to an absolute magnitude 
^break ~ —19.6. A transition at about the same magnitude in 
the galaxy properties has been evidenced by Kauffmann et al. 
(2003) and Baldry et al. (2004). 

An additional interesting outcome of our analysis is that the 
Mstar-M/, relation is already established at redshift z w 1.7, in 
line with the theoretical expectation of the anti-hierarchical 
baryon collapse scenario (Granato et al. 2004). 

Applying the same technique to the local velocity disper- 
sion function of galaxies and to the black hole mass function, 
we have also computed the er-M/, and Mru-M/, relationships. 

The former is quite close to a single power law a oc M^. The 
latter is again a double power law breaking approximately at 
Mh,bieak, corresponding to Mbh ~9x 10 6 M q . The associated 
velocity dispersion, a ~ 88 km s -1 , is very close to the first 
estimate of the critical velocity dispersion for the gas removal 
by SN explosions given by Dekel & Silk (1986), who found 
a critical halo velocity V c m ~ 120 km s" 1 , corresponding to a 
critical velocity dispersion er cr i t ~ 80kms _I . 

As a test of our results, we combined the Mbh-M/, relation 
[eq. ( I16H with the a- Mh relation [eq. J15I 1 : as shown by the 
lower panel of Fig.|6] the resulting Mbh-ct relation is consis- 
tent with the observational data. 

The relationships we have obtained are model-independent 
and can be interpreted in terms of feedback effects by super- 
novae and AGNs in galactic structures. We also presented 
a simple feedback model, which nicely reproduces the ap- 
proximate proportionalities Mbh cx M star oc M/, observed in 
the high mass range, and the break of these relationships at 
M libreak «3 x 10 H M Q . 

At low masses, the M stal — M/, relation derived here (M stal oc 

Mh' 1 ) is steeper than that yielded by the model (M star oc M^ 3 ), 



and would imply that only a tiny fraction of the baryons ini- 
tially associated with the halo remains within it in the form 
of stars (and we know that the gas does not add much to the 
baryon content of galaxies). On the other hand, if the amount 
of stars formed is so low, for a standard Salpeter IMF the 
energy injected by SN explosions is insufficient to expel the 
residual gas if the baryon fraction is close to the cosmic value. 
Thus, if the slope of the M stal — M/, relation really is as steep as 
its face value suggests, we must conclude that either the initial 
baryon fraction in low-mass galaxies was substantially lower 
than the cosmic value (due, e.g., to a pre-heating of the inter- 
galactic medium hampering the infall of baryons into shallow 
potential wells) or that the SN feedback in these objects was 
substantially stronger than in more massive galaxies. As dis- 
cussed in §|3] however, the uncertainties on the low luminosity 
portion of the LF are large enough to allow for a flatter slope, 
closer to the model prediction and almost sufficient to grant 
the gas removal by SN feedback. 

The errors shown in the figures mostly reflect uncertainties 
in the M/L ratio. We must not forget however, other error 
sources. For example, Fig.[3j? shows that different choices for 
the GHMF yield a systematic difference in the results, that, at 
the high-mass end, become comparable to the scatter consid- 
ered in the same figure. Further uncertainties come from esti- 
mates of the GSMF; these are illustrated, in Figs. [3] 0] |5] and 
by comparisons with results obtained using the GSMF by 
Belletal. (2003). Nevertheless, our approach provides results 
consistent with observations, and have comparable uncertain- 
ties. Moreover, since our approach bypasses any assumption 
on the DM profile, it could provide a valuable tool to discrim- 
inate among the different models of DM mass distribution in 
galaxies. 

Our analysis has shown that the relationships presented 
above bear the imprint of the processes ruling the galaxy for- 
mation and evolution. Models should eventually comply with 
them. 
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APPENDIX 
A SIMPLE FEEDBACK MODEL 
The rate at which the gas mass M; n f falls toward the central star forming regions can be written as 

M inf (0 = -^. (Al) 

tc 

The infalling gas mass thus declines exponentially M; n f(f) = M; n f(0) exp(-f /t c ), where t c = max[f coo i(r v i r ),f 1 iyn('"vir)] is the maximum 
between the cooling and the dynamical time at the virial radius, while Mi„f(0) = / cosra M/, is the initial gas mass. 
The time derivative of the cold, star forming gas mass is given by 

M cold (f ) = - m +R^(t) - am , (A2) 

where ip(t) = M stal is the Star Formation Rate (SFR), R is the fraction of mass restituted by evolved stars (R w 0.3 for a Salpeter 
IMF), and 

Nsn^snEsn , 

a = - (A3) 

E B 

is effective efficiency for the removal of cold gas by the supernova feedback. Here Nsn is the number of SNe per unit solar mass 
of condensed stars, esnEsn the energy per SN used to remove the cold gas, and Eg the binding energy of the gas within the DM 
halo, per unit gas mass. Following Zhao et al. (2003) and Mo & Mao (2004), the latter quantity can be written as 

^ = ^V v 2 ir /(c)(l+/ cosm ), (A4) 

where V v i, is the circular velocity at the virial radius for a halo mass M/„ /(c) « 1 is a weak function of the concentration c, and 
we have assumed that, initially, the gas fraction is equal to the cosmic baryon to dark-matter mass density ratio, / cosm . Taking 
into account the dependence of V v i, on the halo mass and redshift, we get, for z > 1: 

E B ^32xl0^( l 4^)(-^L-] (cm/s) 2. (A5) 



4 J V W 2 Mr. 



The effective efficiency is well approximated by 



N SN \ (e_M\ ( E SN \ ( M h ■ 1 : ' 



Further setting 



l2 [s^){u){To^- g ){—) [w^J ' (A6) 

m=^, (a?) 



star 



f star being the star-formation ti mesc ale, eq. iA2\ is easily solved forM co id(0- The mass M stal - cycled through stars is then straight- 
forwardly obtained, using eq. (IA7> : 



AA 1; ,(/)= / ^Vf'=^P 



1 - exp(-f/f c )+ - exp(-i 7 f /t c ) 

57 — 1 57 — 1 



(A8) 



where 7 = 1 -R + a. In the above formula, s = t c /t st . dr ^> 1, since we expect that in the central, clumpy regions the cooling and 
dynamical times are shorter than f c , which is estimated at the virial radius. The mass in stars at the present time only includes the 
fraction / surv of stars still surviving: M" t °™ = / sm v^star(fnow)- The survived fraction depends on the IMF and on the history of star 
formation; as a reference, for a Salpeter IMF after about 10 Gyr from a burst we have / surv w 0.6. If we assume that most of the 
stellar feedback comes from SN explosions, then the total energy injected into the gas is given by 

17 e 1 at ha o m58 f e SN\ ( E SN \ ( N SN \ M stal - 

£star = ^snEsnNsnM.,.,,- « 8 x 10 — — r — erg. (A9) 

star SN SN SN stai V 0. 1 / \ 10 er g / V 8 X 10 / 10 M 

As long as a significant amount of cool gas is present in the central regions, we can imagine that there are mechanisms able to 
remove angular momentum from gas clouds bringing them into a reservoir around the central BH. One of these mechanisms is 
the radiation drag (Kawakatu & Umemura 2002) according to which, as shown by Granato et al. (2004), the reservoir is fuelled 
at a rate 

M res = 1.2 x 10"V(0(1 ~e~ T ) , (A10) 

where r is the effective optical depth of the central star forming regions [cf. Eqs. (14), (15) and (16) of Granato et al. 2004]. If 
most of the mass in the reservoir is ultimately accreted onto the central BH, we expect 

M BH «1.2xl0- 3 M star (l-O- (All) 
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1 /3 

Granato et al. (2004) assumed that the effective optical depth depends on the cold gas metallicity and mass r oc ZM g a S . The 
outcome of their numerical code yields, on average, ZaMjJ 3 in the mass range 10 n M Q < M h < 3 x 10 13 M Q [cf. their Figs. (5) 
and (8)]. Since M gas ~ f cosm M h , one gets 

TocM 2 h /3 . (A12) 

As for the AGN feedback, we use the prescription of Granato et al. (2004); we rewrite their eq. (28) for the kinetic luminosity 
that can be extracted from AGN-driven winds, as 

L K = » 4.2 x 10 44 / C JV22 ergs" 1 , (A13) 

where f c is the covering factor of the wind and N22 is the gas column density in units of 10 22 cm" 2 . If we assume that the BH mass 
is growing at around the Eddington rate, the total kinetic energy in winds emitted by a BH of mass Mbh is Ek ~ (2 /3) ^/^(Mbh)- 
This shows that the action of the AGN occurs on a short timescale, the e-folding timescale t e f = (e?£)/(l - e) (where ?£ is the 
Eddington time and e is the BH mass to energy conversion efficiency; for e = 0.U e /~5x 10 7 yr). 
If a fraction /}, of the AGN kinetic energy is transferred to the gas, its total energy input is 

^G N = /^^3.6xlO-/ ft/c AT 22 ( T l-)(^^) (^)V (A14) 

Since studies of BAL QSOs suggest that N22 > 30, f c > 0.1 (see, e.g., Chartas et al. 2002; Chartas, Brandt & Gallagher 2003) 
and fh > 0.3 (see, e.g., Inoue & Sasaki 2001; Nath & Roychowdhury 2002), we can take fhfcNn w 1. It is interesting to compare 
this energy input to the total energy released by accretion, E acc = 1.8 x 10 62 (e/l -£)(M B h/10 8 Mq) erg: 

-2, (Jinn 



E AGN w 2 x IQ-'E, CC I } (A15) 



